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SUMMARY 

. A solution for the compressible potential flow past slender bodies 
of revolution has been derived by an iteration procedure similar to 
that of the Rayleigh-Janzen and Prandtl-Ackeret methods . The solution 
has been analyzed with respect to transonic similarity. The results 
obtained are in approximate agreement with those of von Karman in the 
region of the flow field not too close to the body. In the neighbor- 
hood of the body, a different similarity law is obtained. This new 
similarity law holds for variations in thickness ratio and Mach number, 
but not for variations in specific-heat ratio. In addition, this law 
appears to be limited in applicability to extremely slender bodies of 
revolution probably outside the range of practical interest. The dif- 
ferences between the results of the present investigation and those of 
von Karman are interpreted in terms of the manner in which the boundary 
condition on the body is satisfied and of the nature of the singularity 
of the solution near the axis. 


INTRODUCTION 

Transonic similarity rules for thin airfoils and slender bodies of 
revolution have been derived by von Karman (reference l) . These rules 
for the case of two-dimensional flow are verified in reference 2 by an 
iteration procedure similar to that of the Rayleigh-Janzen method (ref- 
erence 3) and the Prandtl-Ackeret method (reference 4). 

The analogous investigation for bodies of revolution was made at 
the NACA Lewis laboratory and is presented herein. A solution for the 
compressible potential flow past a slender body of revolution is obtained 
by the same iteration and transonic limiting procedure used in reference 2. 
As in reference 2, in each step of the iteration procedure the boundary 
conditions were satisfied on the body. 

The solution obtained yielded transonic similarity rules that are,, 
in part, different from those of reference 1. The differences appear 
to result from the. method of satisfying the boundary condition on the body 
and from the nature of the singularity of the solution near the axis. 
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As in the two-dimensional investigation, the influence of stag- 
nation points has not been considered herein. 


ANALYSIS 

General formulation. - The partial differential equation for a 
compressible, axially symmetric, isentropic and irrotational flow with 
free-stream velocity U in cylindrical coordinates x,r ,0 (fig. 1) 
is 


[a 2 - (U+u) 2 ]^ + (a 2 -v 2 )cp rr - 2(U+u)cp r cp xr + a 2 ^ =0 (1) 

in which the following notation has been used: 
a local speed of sound 
U+u resultant velocity in x-direction 
v resultant velocity in r-direction 

cp perturbation velocity potential defined by u = cp^, y - cp^ 

(Subscripts denote differentiation with respect to variable noted.) 

(All symbols used herein are defined in the appendix.) The local speed 
of sound a is related to the free-stream speed of sound a , the 

ratio of specific heats y, and the local velocity by the Bernoulli 
equation 

a 2 . a Q 2 - Zzi (2Uu + u 2 + v 2 ) (2) 


In accordance with the Prandtl-Ackert type of procedure, 
equation (l) will be written in a form in which the linear terms appear 
on the left side of the equation and the nonlinear terms on the right 
side. A solution will then be sought in the range of free-stream Mach 
number close to 1 and thickness ratio close to 0 on the assumption that 
the flow pattern . obtained by inclusion of the nonlinear terms will differ 
by only a small, amount from that obtained with only the linear terms. 

This assumption will then be made plausible by the resulting form of the 
solution. The coefficient of cp^ in equation (1). is therefore expressed 
with the aid of equation (2), in the form 


a* - (U + u ) 2 . - "o 2 K 0 <(ffj 


( 3 ) 
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where 

Mq free-stream Mach number, u/eq 

3 2 = 1-M 0 2 (4a) 

r M - «0 2 ( x + M o 2 ) <«>> 


For convenience, the free-stream velocity is taken as the unit velocity 
so that u/u, v/u, and cp/u may be written as u, v, and cp, 
respectively. The differential equation (l) can now be expressed in the 
form 


- [ r M< a f , x + 'f I 2 ) + ^ “oV 2 ] 'Pet + (5) 


The boundary conditions of the problem are: (1) The perturbation 
velocities vanish at infinity, and (2) the flow follows the contour of 
the body. Thus, at infinity. 


= <P r = 0 (6) 


and on the body, 

- Tg(x) (7a) 

cp r = T(l4cp x ) ^(x) (7b) 

where 

T lateral distance ratio of body 

g(x) function characterizing shape of body and of order of magnitude 1 
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All lengths are expressed in terms of the chord of the body as 1. 

In order- to obtain the Laplacian of cp on the left side of 
equation (5), the affine transformation 

0) = pr 

is introduced. The differential equation (5) becomes 

Aqp « ^ ^, 2 (2Cp 3c -rtp JC 2 +p\ 0 2 )Aq) + 



|pf (2Cp x +qp x 2) + ^ V%) 2 J ^xx + M o 2 P 2c P(0 2< Pa)(C + 

2M 0 2 ( I -Kp 3c )c^ 0 cp !c(0 

where the Laplacian of cp in cylindrical coordinates is 

Ac P = ^xx + ^(00) + HT 

The boundary conditions (7a) and (7b) become, 
at infinity, - 

*X - %> - 0 . 


on the body, . 

V 

“b = T P S(x) 

<P(0- cp x = I ^(x) 


(9) 

( 10 ) 

( 11 ) 

(12a) 

(12b) 


The formulation of the problem so far is exact. A solution is now 
sought that is applicable in the range T and (3 close to zero (t~ 0, 

P ~0). This solution is referred to hereinafter as "the small-perturbation 
transonic limiting solution." In the range of . T and p under considera- 
tion (t ~0, p ~0), the perturbation velocities will be assumed small 
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compared with free -stream velocity or |u|« 1, |v|« 1. Thus, as is 
usual in this type of procedure, the right side of equation (9) is 
considered to produce a small perturbation from the linear case and a 
solution of the system of equations (9), (11), and (12) will be sought 
in the form 


12 3 

<p = cp + cp + cp + 


(13) 


in which each term is of a lesser order of magnitude than the preceding 

12 

one. The following boundary conditions on qj, cp , . . . will accord- 
ingly be taken as the equivalent of the boundary conditions (11) and (12) 

At infinity, 

^Px = < ^(o = ^ n = ’ • * (14) 

On the body, 

% = P T g(x) (15a) 

= if (15b) 

?co - ^ g^x) $ x = °; n/l (15c) 

In order to obtain a solution of the system of equations (9), 

(13), (14), and (15), equation (13) is inserted into the differential 
equation (9) and a typical Laplacian term on the left, such as Acip, is 
equated to the sum of those terms on the right that contain the super- 
script n-1 and that may also contain any of the superscripts n-2, 

n-3, ... 1. The right side of the differential equation for cp con- 
sists of a sum of terms 'of which, for a range of Mq near 1 and j 

near zero, seme will be of highest order of magnitude. The solution 
corresponding to these terms constitutes the small-perturbation tran- 
sonic limiting solution for. $. 

Transformation to prolate-elliptic coordinates . - For the problem 
of the flow past an isolated body of revolution, it is convenient in 
satisfying the boundary conditions to use a system of prolate -elliptic 
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coordinates (reference 5). The transformation from cylindrical 
coordinates x,(0 to elliptic coordinates n,X is given by 


x «'*iX (16a) 

W = -/(X 2 -l)(l-^ 2 ) (16b) 


The surfaces X = constant, p. = constant are confocal ellipsoids and 
hyperboloids of two sheets, respectively, the common foci being at 
x = +1, (0=0. The values of X may range from 1 to infinity, whereas p 
varies between -1 and +1. 

The Laplacian of cp is, in elliptic coordinates, . 



Transformations from derivatives with respect to x and (0 to 
derivatives with respect to p and X will be needed for subsequent 
analysis. From the transformation equation (16), these relations are 


^x - 


n(* 2 -l) 

(X 2 -p 2 ) 




X(l-u 2 ) 

(X 2 -n 2 ) 


% 


(18) 


Cp. _ ^ 2 -D(i-F -!l 

^(0 - 9 .\ 


(X 2 -H 2 ) 


(^X-^m) 


(19) 


(* -u ) L J • (xV) 


2kX(1-H 2 )(X 2 -1) t X 2 (1-h 2 ) 2 

,.2 2.2 .V i.'S "W 

(A ) 


( x 2 -,, 2 ) 1 


( 20 ) 
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<P« 


COO) - 


(\ 2 2\ 3 
(A -|i ) 


( X 2 -2n 2 X 2 -Hi 2 ) ( 1 -|i 2 ) «ii 2 (X 2 -1 ) ( 2 -X 2 


K 


[( X 2 t2h 2 X 2 -^ 2 ) ( X 2 -1 ) -X 2 ( 1-kt 2 ) ( 2 -X 2 -)i 2 )] MCp^ 


> + 


(X 2 -] 




(X 2 -,,") 


' (21) 




XOD 


if^-Dd-n 2 ) 

r M 

(X 2 (3-X 2 ) + |^ 2 (1-3X 2 ) 

<*v> 2 

L 

(X 2 -m 2 ) J 


<PX - 


U 2 (5-u 2 ) 


n 2 ) +xW)~ | 
(x% 2 ) J 11 


(iX(X 2 -l)cp-^^ + (XW^V)^ - nX(l-n 2 )cp^l 


( 22 ) 


The variable X has the limiting value 1 on a body of revolution 
as the thickness ratio T of the body approaches zero. It will there- 
fore be convenient to have the limiting forms , for X -*-1, of 
equations (18) to (22). Defining a new variable^ l by 


l = X 2 -l (23) 

and denoting by ~ the limiting form for X-*> 1, equations (18) to 
(22) take on the limiting forms 
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z i/2 

9(0 71 75 (2 V^V 

(i-n 2 ) 


(25) 


^xx ~ 2 2 ["2(1+34 ) “ ^( 3+4 ) Zcp^ + 44 Z 2 

(1-4 ) L 


^0X0 ~ 


dV) 


44(l-4 2 )Zcp^ + (1“4 2 ) 2 CP J 
2CP Z - 4CP U + l (4C^ Z - 44cp^ + 4 2 cp 44 )j 


(26) 


(27) 


cp. 


* i 1 ^ 2 r 2 

~ ~ 2 3/2 L^z " (1+ ^ + 4 ^zz + 2(1 ^' 


XO) 


(1-4.^) 


)c Plu - . 

(28) 


The boundary condition (14) becomes, in elliptic coordinates, at 
infinity, 


(x 2 -* 2 ) (x 2 -h 2 ) 


| n = 1,2, 3, . . . 


V i X 2 .l) (1 V) ' ^ n , _ 0 
(A 2 -4 2 ) A ^ 


(29) 


The boundary condition (15) has, as a limiting form in elliptic coordi- 
nates, on the body. 


Z 


b " 


Mi. 

(l-M 2 ) 


(30a) 
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-TJTJZ <2 ^^ > - f «l (:t) (^2 - P ^ (x) (30b) 

l 1 / 2 n n , x . , /2pl n n\ . . . 

P72 ( 2 V%) " p g 3c (x) Vl+%) = 0; n / 1 (30c) 


First approximation . - The solution for the first approximation <p 
satisfies the homogeneous Laplace equation 

[(l-p 2 )^ + [(X 2 -!)^- 0 (31) 


subject to the boundary conditions (29), (30a), and (30b). A solution 
of equation (31) appropriate to the present problem is given by (ref- ’ 
erence 5) 


^i(u', X) 



L QnU) P>) 


(32) 


vhere P n (p) and (^(A) are the Legendre functions of the first and 
second kind, respectively, and the An are constants . The function 
P n (|i) is a polynomial in p of degree n, free of singularities for 
finite values of p and is given by (reference 6) 


P 0 (b) - 1 
= M 

P 2 (p) = | (3p 2 -l) 


P>) 


1 d p , = ^ (-l) J (2n-2j)'. n -2J 

2 n ni dp n * ~ pg' 2 n j'.(n-j)'.(n-2j)l 


(33) 


I 
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where m « n/2 for n even and m = .(n-l)/2 for n v odd. The 
function <^(X) has a logarithmic singularity at X = ±1 and is 
given "by (reference 7) 


Qq(X) = i log 


1 (X+l) 2 

X 2 -l 


Q X (X) = -4^- los - 1 


X 2 -l 


^(X) = |P 2 (X) log - f P X (X) 


X 2 -l 2 
dX 


00 

= P n (M r “ g r =T2 

\ (X*-ijLp b (mJ 

- I F n< X > ** ^ ? n-l (X) - fffer P - 3(X) 


(34) 

For large values of X, the function ^(X) has, as limiting form, 
(reference 6) 




2 n (n’) 2 
(2n+l) lX n+1 


The function Q, n ( X) and its derivatives therefore approach zero as X 
"becomes infinite, so that the boundary condition at infinity 
(equation (29)) is satisfied. The constants 3^ in equation (32) are 

determined so that the boundary condition at the body (equations (30a) 
and (30b)) is satisfied. In the neighborhood of slender bodies, X 
approaches 1 so that the limiting forms of Q^X) and its derivatives 
are applicable near the body. These are, from equations (34), 

% 

Qjl) ~ - \ log l (35a) 
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(35b) 


~ (35c) 

2Z 2 

vhere the primes are used to denote differentiation with respect to l. 
Primes will be used, where convenient, to denote differentiation when 
only one variable is Involved. The limiting form of the function 0^(1) 
is evidently Independent of the subscript n, so that equation (32) has 
a limiting form In the neighborhood of the body 

•Wj I) = - \ log \ P n (n) ' (36) 

n^O 

Insertion of equation (36) Into the body boundary condition (30b) for 

1 

the determination of ^ yields the equation 



For convenience in discussion, identifying supernumerals have been placed 
above the various terms in equation (37). These terms will next be 
compared in the small-perturbation transonic limit. The ratio of term 
<D to term (I) is, by equation (30a), of order T^p^ log TP, whereas 
the ratio of term © to term @ is of order l/logTP. Both of these 
ratios have a limiting value of zero regardless of how t and p 
individually approach zero. It is therefore permissible to neglect © 
with respect to © and ( 5 ) with respect to @ . The ratio of term (|) 
to term (p is of order log Tp; the limiting value of this ratio 
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depends, on the manner in which T and p approach zero. Although 
T 2 log Tp will, in general, he small for 3mall t, values of p 
exist that, for small but fixed values of' T, will make t 2 log tP 
comparable to 1. The quantities T and p are to be regarded as 
independent variables in going to the small -perturbation transonic 
limit; i,t is therefore necessary to retain term (4) along with term 
(l) in equation (37). Term (4) arises from the Cf^. term in the boundary 
condition (12) and is generally neglected in small -perturbation analysis. 

When terms (2) and (3) are neglected, and log ^/(l-n 2 ) 1 / 2 is neglected 
relative to log tp, equation (37) reduces to 


v p >> 




Kgs x> 


2 

‘ T 8«x 


(38) 


where 


K = T 2 log TP 


(39) 


Both sides of equation (38) may be expressed as series of ortho- 
gonal Legendre polynomials in n and the coefficients of corresponding 
Legendre polynomials equated to each other. The result is a set of 

simultaneous linear equations for the A. , the solution of which may be 
indicated as 


1 


T 2 ^Osb) 


(40) 


where a^Kjb) denotes a function of the parameter K and of body 
shape; the functions g. and g^ are of order of magnitude 1 so 

that it is plausible to assume is also of order 1. 

The solution for the first approximation $ is then 

b = t 2 ^ » n Q^X) P n (|i) (41) 

n=0 

Near the body qj has the form (from equation (35a)) 
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1 

cp 


log_Z 


00 „ 


®a 


n=0 


(42) 


The perturbation velocity cp x near the body is, from equations (24) 
and (42), • 


1- 



1 

a 

n 


P n (n) + \ log l 



(43) 


On the body or in the flow field not too far from the body, the second 
term in the bracket dominates 'the first term in the transonic range 
P ~ 0. In particular, on the body cp^ has the form 

» 00 

<p x ~ - T 2 log Tp ^ a n P n '(n) = - it ^ an P n '(n) (44) 

• n=0 h^O 


The radial velocity component A near the body is, from equations (8) 
and (25), 



T hi 1 ' 2 

2^72 


dV) 




The second term in the bracket is negligible compared with the first 
term, so that on the body equation (45) assumes the form 

v r ~ - V 1 *' <«) 

n^O 


In' the flow field, which will be defined as the region for which l 
and log 2 are considered to be of order 1, the distrubance velocity 

qig is obtained from equations (18) and (41), whereas the radial velocity 

q> r is given by equations (8), (19), and (41). Because only the order 

of magnitude of the solution in the flow field, rather than its specific 
form, is desired, the process may be considerably simplified. The 
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functions Q n ( A) and F n (n ) will contribute only terms of order of 

magnitude 1 to the solution. The form of the solution in the flow 
field may therefore be obtained from the solution near the body by 
<■ taking I and log Z of order 1. The solution for cp x is therefore, 
from equation (43), of the form 


<P X ~ t 2 g 1 (K,b) 


(47) 


where g 1 (K,b) is a function of K and of body shape, and is of order 
of magnitude 1. The radial velocity <p ' in the flow field has the foraj 
from equation (45), 


hr ~ ^(^b) 


(48) 


where h^(K,b) is a function of K and of body shape, and is of order 
of magnitude 1. 

1 

The resultant velocity increment A is defined by 


1 

A 




<i > 2 J 


-I 1/2 


- 1 


(49) 


The ratio (cp )^/p is of order T 2 p 2 in the flow field and of order 
r x 1 2 

l/log Tp on the body so that (cp r ) may be neglected. The' resultant 

1 1 
velocity increment A is therefore approximated by 9 X « 

2 

Second approximation . - The second approximation cp is a solution 
of the non-homogeneous Laplace equation 


m q 2 ( 2c P x A 2 +3 2 %) 2 )a<P +■ 


© 


P 2 


© ® 

< 2 V* x 2 > * “o 4 ^ 2 

©© 


1 

9 + 

XX 


M 0 2p Wo) + 2M 0 2 (lA)^XO) 


(50) 
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which satisfies the boundary conditions (29), (30a), and (30c). In 
solving equation (50), only those terms on the right side that are of 
highest order of magnitude in the small -perturbation transonic limit 
t*“*0, 3 — >0 will be retained. The plausible assumption is made in 

this connection that the highest-order part of the solution of a 
complete differential equation is the same as the complete solution of 
the differential equation in which only highest-order terms have been 
retained. For further simplification, an explicit solution will be 
sought only near and on the body, so that the limiting forms for the 
various derivatives in equations (24) to (28) may be used. The sub- 
sequent analysis will therefore only indicate implicitly the form of 
the solution in the flow field. 

The highest-order terms in qj • and <p. are, from equations (43) 
and (45), 

~ - T 2 21 log Z 
x 2 

1 T 2 R 

z 1 / 2 (1- n 2 ) 1 / 2 


where R = R(k ,b,n) = y a^( K,b) P n (n) and the primes denote differ- 

n=0 

entiation with respect to n . The highest-order contributions to the 
second-order derivatives of cp are, by equations (41) and (24) to (28), 


1 2 t)M 

^xx - " T \ 108 1 

1 ^ T 2 R' 

^0) - ^172 (1 ^ Z) l/2 


1 



R 

(l-u 2 ) 


The various terms on the right side of equation (50) will now be 
examined with regard to their orders of magnitude. Term (l) vanishes 
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by virtue of the solution for the first approximation. The ratio of 

term @ to term ( 2 ) is of order t 2 3 2 / Z log l and approaches zero 
both in the flow, field and on the body. It is therefore permissible to 
neglect term ( 3 ) . The ratio of term (2) or ( 5 ) to term, ( 3 ) is of order 

0 2 /Z log 7, which is negligible in the flow field; near the body, the 
ratio may be small or large depending on the manner in which T — » 0 and 

p— >0. Terms (2) and ( 5 ) must therefore be retained and equation (50) 

2 

for the second approximation cp reduces to 


^ p 

Acp i M 

t 4 V 


R'R" log 2 Z - g- (B' ) 2 R"T 2 log 3 Z^ + 


Mo 2 B 3 t 2 p 2 


dV ) 2 * 2 


2Mrt 2 - RR ' i. . m2 BCR') 2 T 2 log l 
2M ° (W«) » * (l-, 2 ) * 


(51) 


A typical term on the right side of equation (51) is, aside from 
factors depending on the physical parameters T, 7 , and Mq, the 

limiting form for A — *1 of some function of p. and A, which may be 
denoted by f(n,A). In order to solve equation (51), it is assumed that 

A O 

the product of f(u,A) and (A c -n ) may be expanded into a series of 
orthogonal Legendre polynomials. The right side of equation (51) may 
contain singularities at the end points p * +1 because of the factor 

— 1 — , which becomes infinite at these points. It may therefore be 
1-H 2 

necessary to exclude small but finite neighborhoods about n e ±1 from 
the interval of expansion and consequently from the domain of definition 
2 

of cp. At this point the assumption is made that the effect of the stag* 
nation points, if they occur, is negligible. The expansion for f(p.,A) 
may then be expressed by 


(X 2 -u 2 ) f(n,») = 



% V*) p n ^ ) 


(52) 


where depends on body shape and h n is a function of A. Then the 

solution of equation (51) corresponding to the term f(n,A) is, from 
equation (17), the solution of the non-homogeneous equation 
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(x 2 -n 2 )Acp = 




Si V X) P n^>. ( 53 ) 


A solution of equation (53) is now sought in the form 


q> = 



(54) 


where q Q (X) is a function of X to he determined. 

Inserting equation (54) in equation (53) and noting that P n (n) 
satisfies Legendre's differential equation yields, when coefficients of 
P n (|i) are equated, the non-homogeneous Legendre differential equation 

[(X 2 -l) q Q '(X)] r n(n+l) q n (X) = h^X) (55) 

The solution of equation (55) is the sum of the complementary solution 
taken proportional to Qq(X) to satisfy the boundary condition at 

infinity ancLthe particular integral r n (X), which may he expressed in 
terms of two indefinite integrals hy 

, * <M 

The function r n satisfies the houndary conditions at infinity. The 
solution of equation (53) may therefore he written as 


hn(X) 


P n (X) d\ 


d\ 


(56) 


Cp s 



^ p n^) + 



% r n< X ) p n^) 


( 57 ) 


The limiting form of equation (56) for X — ►! or l — >0 is 
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The subscript n is omitted from r(Z) inasmuch as h Q (A) is 

independent of the subscript n in the limit Z — >0. Equation (57) 
then becomes, for I — >0, 

ep~ T 4 - i iog An P n (*x) + r(l) ^ P n (n) (59) 

L n=0 n=0 J 

The quantity r(Z) may now be evaluated for each term on the right side 
of equation (51) and the solution for cp Consequently is 


log Z^ ^ P n (^)+ | (1-p 2 ) -*rR’B”Z log 2 Z- ^{R' ) 2 R"t 2 Z log 3 z] + 
T ru=0 ^ 


i*o 2 


+ iM. 

(l-n 2 ) Z 4 o 


HR' log^ Z- 


1_ 

24 


Mo 2 R(R‘) 2 t 2 


log 5 Z (60) 


Equation (60) explicitly indicates the solution only near the body, 
where Z ~0. In the flow field (Z ~ 1 and log Z ~ 1), equation (60) 
indicates the solution only as to order of magnitude, as previously 
mentioned. This limitation, however, is unimportant in the investigation. 
The limiting form of the solution as T ~ 0 (for example, equation (60)) 
is used herein because of the resultant simplification in the analysis. 

The constants ^ in equation (60) are determined from the boundary con- 
ditions (30a) and (30c). Insertion of equation (60) into equation (30c) 
yields • 
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2 R 2 2 

where Z ia to he evaluated at ^ 8 

l~p . 2 

equation (61) become negligible compared with the checked terms in the 


The unchecked terms in 


small -perturbat ion transonic limit. Equation (61) therefore reduces to 


^ Xi P n (n) - kggx ^ P n !‘(n) 




= log Tp f2M 0 2 ER' + K ^r M g 2 R'E"-M 0 2 R(E') 2 -M 0 2 gg x (RB')’] - 


g 2 (R’) 2 R" - |m 0 2 ggx 


[«.■,§) 


(62) 


Equation (62) may be solved for the coefficients A- in the same manner 

1 

as was indicated for the coefficients A n occurring in the first 

approximation <p. As in equation (52)/ it may be necessary to exclude 
the region around the- end points p = +1. The result may be expressed 
as 


2 2 

A^ = log Tp a n (K,MQ J Fj ( j,b) 


where aj/^Mg^^b) denotes a function of the parameters K, ' Mg, . 
and the body shape and is of the order of magnitude 1. 


Although the term in IT cannot be neglected in equation (51), it. 


can be neglected in equation (60) and in the boundary condition 


(equation (61)). Near the body, the solution for cp given by 


equation (60) therefore becomes, when the R° term is neglected, 

2 CO 2 r 

~ 1 log TP log z£ Xx P n (^) +|(1-^ 2 )-^[r , R ,, Z log 2 Z -^(R’) 2 R"T 2 log 3 z] + 


n=0 


y Mq 2 !®- log 2 Z - ^ Mg 2 R(E') 2 T 2 log 3 Z 


(63) 
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An explicit solution for the potential in the flow field cannot he 
obtained from equation (63) but, as previously mentioned, the order of 
magnitude of the solution in this region may be inferred from 
equation (63) by taking l and log l of order 1. In the flow field, 

2 2 

the disturbance velocities cp. x and cpj. thus have the form 


I, ~ T 2 (c-LK+CgC) 


(64) 


<P r = ~ T2 P( d i K+d 2 € ) 


(65) 


where e = 


T 2 r, 


p 


c 1 • and d 1 are functions of K, Mq, T m , and body 


shape, and C 2 and d 2 are functions of only K and body shape. The 
c's and d's . are of order of magnitude 1. 

In the neighborhood of the body (equation (30a)), the quantity l 
begins to contribute to the order of magnitude of the terms in 
equation (63). The disturbance velocities on the body are thus given 
by 


d. p 

cp x ~ K K 


x vM + - 1 % 2 [r(k') 2 ]k 


u *rr 

Cp r ~ - - K< 
r g 



(66) 


X. 


2M q 2 RR' 


[mo 2 R(R’ 


)2 -r M g 2 R’R^ 


K + 


ir^R 1 ) 2 !*" K 


(67) 


Aside from the dependence of the constants a on K , the dominant 
2 n ’ 

terms contributing to cp x on the body (equation (66)), come from the 

%&a> and terms in the differential equation (50). The first 

w 300 2 11 
term of cp x in the flow field (equation (64)) results from the cp m_ ( . 

term. The second term of equation (64) comes from the C P X C P 0; term in 
(2) in equation (50). The expressions for cpp on the body and in the 

flow field (equation (67) and (65), respectively) may be similarly 
^ analyzed. 
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DISCUSSION 

Transonic similarity . - The results obtained thus far in the first 
two approximations may be summarized as follows : 


On the body, 


- K 

R + (S-M 0 2 RR')* + | M q 2 R(R' ) 2 

(68) 

- 

Jk' +[^S'-M 0 2 (EE')^K + |m q 2 £r(E') 2 ]' Kj 

► (69) 

,xJ 

g 

R+(S-2Mq 2 ER' )k + |m 0 2 R(R') 2 - r M g 2 R'R""J 

K 2 + ir M g 2 (R') 2 R"K 3 , 


- 

(70) 


where 

oo 

® ~ ^ ' ®n ^n^^ 
n=0 

2 

S=^a n P n („) 

n=0 

In the flow field (Z ~ 1, log Z - l), 

cp ~ T 2 (b x + b 2 K+b 3 £ ) (71) 

<VT 2 (c 1+ c 2 k + c 3 € ) (72) 

cp r ~ t 2 3 (d 1 +d 2 K^ 3 £ ) 


(73) 
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where bg, Cg, and dg are functions of K, Mq, T m , and body shape 

and the other b's, c's, and d's are functions of only K and body 
shape. The b's, c's, and d's are of order of magnitude 1. It 
appears plausible to assume that higher approximations would not alter 
the results obtained thus far concerning the dependence of the potential 
on the parameters K and €. 

The solution given by equations (68) to (73) will next be considered 
from the viewpoint of transonic similarity; that is, the possible 
dependence of the solution on less than three combinations of the 
physical parameters T, Mg, and y in the small-perturbation tran- 
sonic limit T — > 0 , 6 — > 0 will be investigated. In the limit T — >0, 

8 — K), the solution for the body (equations (68) to (70)) becomes a 
function of the parameter k, the ratio of specific heats y, and the 
body shape. Hence, for constant y, a similarity rule exists on the 
body with respect to variations in t and 6 through the similarity 
parameter K. In the derivation of this result, it was necessary to 
neglect terms of order 1 in comparison with log l; that is, in compari- 
son with log t| 3 on the body. The quantity l , however, must be 
extremely small before log l begins to dominate terms of order 1. The 
foregoing similarity rule may therefore possibly be limited to extremely 
slender bodies with thickness ratios not in the range of practical 
interest. If the body is not extremely slender, the potential in the 
neighborhood of the body depends not only on the parameter K but also 
on the thickness ratio T in a complicated manner. (See, for example, 
equations (41) and (63).) 

For the flow-field solution, (equations (71) to (73)), in the limit 
T — * 0, p — >0, the coefficients bg, c 2 , and dg become functions of 
K, 7, and body shape, whereas the other b's, c's, and d's remain 
functions of k and body shape. The flow-field solution thus becomes, 
for T — »0, p — >0, a function of three parameters k, c, and y, so 
that apparently no simple similarity law exists in the flow field. The 
original physical parameters of the problem (t, M q , and y) have 
not been reduced in number, so that no apparent simplification of the 
dependence of the potential upon the physical parameters has been 
achieved. Figure 2, which presents curves of |k| /e against 8 for 
several values of the thickness ratio T, shows, however, that the 
parameter £ is much larger than the parameter K for small values of 
8- It then appears desirable to consider those flow-field solutions in 
which € is of order of magnitude 1 and k is negligible compared with 
£ . 


The essential modifications in the analysis previously given are 
that the term in c p x in the body boundary condition (7b) may now be 
neglected in view of equation (38) and that only the terms in C P X C P XX > 
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qyfcu and cp^ qf^ w need be considered on the right side of the differ- 
ential equation (9). The term in c^cp^. is the dominant term in the 

O 

flow field, whereas the c ^o c P x oo an ^ ^0) %)U) ^ erms 8X9 dominant in the 
neighborhood of the body. It is therefore necessary to retain the 
terms c^cp x ^ and for the body boundary condition in order to 

obtain the constants A^ of the complementary solution. 

By carrying through the analysis directly or by inspection of the 
solution equations (71) to ( 73 ), the potential in the flow field may be 
expressed in terms of the single parameter e. Because of the simplifi- 
cations resulting from the condition K « e, this analysis for the 
solution in the flow field has been carried through the first three 
approximations. In the third approximation, as well as in the second 
approximation (equation (63)), the complementary solution becomes 
negligible compared with the particular solution corresponding to the 
term in the flow field. Thus the results for the flow -fie Id 
region to three approximations may be expressed as 


F = fj£ + f 2 € 2 + f 3 e 3 (74) 

where 

r M 

F = -|cp (75) 

and f-p fg, and f 3 are functions only of body shape and are of 
order of magnitude 1. 

In establishing equation (74), it was unnecessary to make the strong 
assumption log Z»1 but rather the much weaker ones, T^«l, (3^«1 

and log Tg ~ 1. It is therefore probable that the similarity law for 
the flow field given by equation (74) holds for a wider range of thickness 
ratios than the similarity law for the body (equation (68)). The 
existence of different similarity laws for the flow field and for the 
body would indicate a transition region where a more complicated relation 
obtains . 

Comparison with reference 1 . - The results obtained in this investi- 
gation are somewhat different from those of reference 1, where a single 
similarity parameter equivalent to the parameter e is derived for both 
the flow field and the body. These differences may perhaps be understood 
as follows : 
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In reference 1, the boundary condition at the body is written in a 
form essentially equivalent to 

rcp r = j 2 SSx tor r >0 
or 


co Fa = egg*. 


for <o— * 0 


(76) 


and nonlinear terms in the differential equation for the potential other 
than the cp y cp YY 'term are neglected. When only the cp^cp^ term on the 

right side is retained and the transformed potential F is used, 
equation (9) becomes 

AF = 2F X F XX ( 77 ) 

If it is assumed that the boundary condition (76) may be evaluated on 
the axis = 0 rather than on the body 0)^ = Tj3 g(x), then the only 

parameter entering the differential equation (77) and the boundary 
condition (76) is e. The potential F(x,0>) should therefore be 
expressible in terms of the single parameter e. On this basis, the 
transonic similarity rules in terms of the parameter c are obtained • 
in' reference 1.. 

The present analysis differs from that of reference 1 in three 
respects. First, the boundary condition (26a), or, more generally, the 
boundary condition 


= T n+1 g n g x (78) 

is, in the present analysis, satisfied on the body, as 

= Tg(x) (79) 

. rather than near the axis , as 

r — >0 (80) 

The use of equations (78) and (79) will evidently yield the same results 
for any value of the exponent n. Equation (80), however, is equivalent 
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to equation (79) only if n = 1, and then only in the first approxi- 
mation; for in the first approximation (equation (41)), cp ~ log r, 
so that repp has a finite nonzero limit as r — *0. Hence, reference 1 

is correct insofar as the first approximation obtained herein is 
concerned. 

In the second approximation (equation (63)), however, the dominant 
2 2 

term is of order log r so that rqp r becomes infinite for r — » 0. 

This result, however, does not invalidate the procedure as regards the 
flow field because the boundary condition (76) is needed only to obtain 
the complementary solution in equation (63). As previously noted, the 
complementary solution is negligible with respect to the particular 
solution corresponding to the F^F^ term in the flow field, at least 
to the first three approximations. 

A second difference between the present analysis and that of 
reference 1 stems from the logarithmic -type singularity that the 
potential cp(x,r) exhibits as r— >0. This singularity affects not 
only the numerical value of the potential cp(x,r) and the . velocity at 
the body r^ =Tg(x), but also its order of magnitude. Hence, even 
after use of equation (79) in satisfying the boundary condition, a 
further use of equation (79) must be made in the solution to determine 
the parameters on which the velocity at the body depends. 

A third difference between the present analysis and that of refer- 
ence 1 is that the singularity of Cp(x,r) as r — » 0 causes the terms 

ln Vxco 811(1 %)^P(00) ( rather than the Cf^cp^ term) in the differ- 
ential equation (9) to be the dominant terms in the neighborhood of the 
body in higher approximations. These terms lead to t}ie similarity para- 
meter K in the present analysis. 

The foregoing analysis indicates that there is no single transonic 
similarity rule for bodies of revolution for both the flow field and 
the body, and that the similarity rule for the body may be limited to 
extremely slender bodies. The use of transonic similarity for bodies 
of revolution may therefore be somewhat limited on a practical basis. 

It is well known, however, that the compressibility effects for a body 
of revolution are much smaller than for the corresponding two-dimensional 
profile (references 8 to 10 ), which may easily be seen fr om 
equations (43) to (45) and (63), (66), and (67) where it is noted that, 
on the -body, the compressibility factor p occurs in the dominant terms 
only through powers of the slowly varying function log Tp . If the 
function log Tp is considered of order 1, then the solution for 
2 

cp is of order T and should be quite small compared with the solution 
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X O 

for cp, which is of order - T . The Prandtl-Glauert rule for todies of 
revolution may therefore he expected to hold for a much wider range of 
free-stream sutsonic Mach numbers them the corresponding rule for two- 
dimensional bodies . 


Lewis Flight Propulsion Laboratory, 

National Advisory Committee for Aeronautics, 
Cleveland, Ohio, August 31, 1950. 
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APPENDIX - SYMBOLS ' . 

The following symbols are used in this report: 
a local speed of sound 

Bq speed of sound in free stream 

r M ' 

F transformed velocity potential, — qp . 

3 2 

g(x) function characterizing shape of body 
Z variable, >\ 2 -l 

Mq free-stream Mach number, u/a^ 

P n Legendre function of first kind » 

On Legendre function of second kind 

R function of K, body shape, and jj. 

r n particular integral 

U free-stream velocity 

u disturbance velocity in x -direction, cp x 

v disturbance velocity in r-direction, cp r 

x,r,6 cylindrical coordinates 

(3 compressibility factor, <\/i-Mq 2 , (equation (4a)) 

r M = M o 2 + ^IT **0^ ’ (nation (4b)) 


7 

A 

e 

K 

A 


ratio of specific heats 
Laplacian 

similarity parameter in flow field, 
similarity parameter in neighborhood 
resultant velocity increment on body 


t2i m 

P 2 

of body, 


T 2 


log Tp 
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prolate -elliptic coordinates for body of revolution 
T lateral distance ratio 

cp perturbation velocity potential 

co transformed r-coordinate, pr 

Subscript: 
b on the body 
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Ratio of absolute value of parameter tc to parameter c, |k | / c 
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Compressibility factor, £ 

Figure 2. - Curves of ratio of absolute value of parameter k to 
parameter e against compressibility factor {3 for several 
thickness ratios. 
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